Sub-Doppler Laser Cooling and Magnetic Trapping of Erbium 
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We investigate cooling mechanisms in magneto-optically and magnetically trapped erbium. We 
find efficient sub-Doppler cooling in our trap, which can persist even in large magnetic fields due to 
the near degeneracy of two Lande g factors. Furthermore, a continuously loaded magnetic trap is 
demonstrated where we observe temperatures below 25 /iK. These favorable cooling and trapping 
properties suggest a number of scientific possibilities for rare-earth atomic physics, including nar- 
row linewidth laser cooling and spectroscopy, unique collision studies, and degenerate bosonic and 
fermionic gases with long-range magnetic dipole coupling. 



I. INTRODUCTION 

The maturation of laser cooling and trapping tech- 
niques has enabled many recent advances with elements 
ranging over much of the periodic table. Narrow inter- 
combination lines in Sr provide ultra-precise atomic fre- 
quency standards [J, B: Yb has also been proposed as 
a frequency standard [3, [3| , with quantum degeneracy 
recently reached in both bosonic [5| and fermionic iso- 
topes [Q]; a Ra magneto-optical trap (MOT) promises 
new tests of fundamental symmetries, while also exhibit- 
ing an interesting blackbody radiation repumping mech- 
anism 0] ; Cr has a large ground-state magnetic moment 
(6/^s)i which is expected to give long-range order to a 
Bose-Einstein condensate [s]. In an extension of laser 
cooling techniques to a highly magnetic rare earth atom, 
one of us recently reported an Er MOT where a novel 
recycling mechanism allows trapping despite significant 
leakage into energy levels that are dark to the cooling 
light jg]. 

Along with the other rare-earth elements, erbium offers 
an exciting combination of favorable properties for future 
cold atomic physics experiments. All six stable isotopes 
have been trapped, with significant populations in three 
bosonic (^^'^Er, ^'^^Er, ^^Er) and one fermionic (^''^Er, 
nuclear spin / = 7/2) species. Er exhibits a variety of 
accessible transitions that are useful for spectroscopy or 
laser cooling, including a broad (36 MHz) line at 401 nm, 
a narrow (8 kHz) line at 841 nm, and an ultra-narrow 
(2 Hz) line in the telecommunication spectrum at 1299 
nm [l3,[ill,[i3|- It has a very large ground-state magnetic 
moment (T/is), which produces a strong magnetic dipole 
interaction at low temperatures. Its closed 5s^ and 6s^ 
electron shells shield the 4/ core electrons, resulting in 
novel scattering properties while also suppressing unde- 
sirable Zeeman relaxation rates |13i] . Additionally, a de- 
terministic ion implantation source based on an Er MOT 
14| could provide a single-ion analog of a solid-state laser 
while dramatically enhancing the prospects for scalab ility 
in recent rare-earth quantum computing proposals p~5| . 
However, none of these avenues can be pursued without 
first understanding the laser cooling and trapping prop- 
erties of this new system. 

The primary result of this paper is an exhibition of 



cooling and trapping properties, both magnetic and op- 
tical, which bolster the growing list of erbium's desir- 
able features for future cold atomic physics. First, we 
show that efficient sub-Doppler cooling occurs in our 
MOT where temperatures reach as low as 100 fiK, nearly 
an order of magnitude below the Doppler temperature 
{Td ~ 860 ^K). Second, we point out that for the 
401 nm transition, the Lande g factors describing the 
Zeeman shift of the ground- and excited-state magnetic 
sublevels are nearly equal. This near degeneracy of Lande 
g should soften the usual constraint that a'^cr^ polariza- 
tion gradient cooling can only occur in negligible mag- 
netic fields and may have interesting applications for 
atomic beam coUimation or atomic waveguiding along 
magnetic field lines. Finally, we show that a magneti- 
cally trapped atom fraction in our MOT is unexpectedly 
cold, with temperatures falling below 25 /iK, lower than 
expectations based on simple thermodynamic considera- 
tions. We believe these magneto-optical properties, cou- 
pled with the unique atomic properties of the rare-earth 
elements, make Er an exciting candidate for a wide range 
of future experiments. 



II. SUB-DOPPLER COOLING IN THE MOT 

Our apparatus was described in [9] . An effusive atomic 
beam is produced in the horizontal (x) direction by 
heating a sample of Er to 1200°C in a crucible with 
a 1 mm aperture. The beam is subsequently deceler- 
ated in a a~ Zeeman slower [16} then loaded into a 
MOT operating on the strong (r/27r = 36 MHz) tran- 
sition at 401 nm. Cooling and trapping light is de- 
rived from a frequency-doubled Ti:sapphire laser. The 
Zeeman slower beam, MOT beams, and a probe beam 
are derived from the same laser, with frequency offsets 
and relative intensities controlled by acousto-optic mod- 
ulators. All laser frequencies are calibrated and locked 
to a fluorescence probe on the atomic beam, with an 
overall uncertainty of ±2.8 MHz j26j. The quadrupole 
magnetic field at the trap position is the sum of fields 
from our MOT coils together with the Zeeman slower 
and compensation coils, with the largest gradient along 
the direction of gravity (z). Typical field gradients are 
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FIG. 1: (Color online) Schematic level diagram of the Er 
MOT and magnetic trap [^. The excited state of the 401 nm 
transition decays into magnetically trapped metastable states 
at a rate Ri « 1700 s~^. Relaxation to the ground state at 
rate R2 and trap loss at rate -Rloss both occur with rates of 
4 s~^ to 5 s~^. -RLoad varies depending on MOT conditions. 

{dB^,dBy,dB^) = (0.20,0.11,-0.33) T/m with ±5% un- 
certainty. All of the measurements described here were 
performed on ^^^Er, with trapped atom populations be- 
tween 10'' and 2 x 10^. 

Er exhibits a complex level structure, with 110 elec- 
tronic states lying between the ground and excited state 
of the 401 nm laser-cooling transition. As shown in Fig.[Tl 
the MOT remains functional because decay of excited 
atoms into metastable intermediate states does not nec- 
essarily result in the loss of those atoms. On the contrary, 
a large fraction of these metastable, dark-state atoms re- 
main magnetically trapped and eventually relax back to 
the ground state. We begin by studying the tempera- 
ture of the MOT, which we measure by observing the 
ballistic expansion of the atom cloud after all optical and 
magnetic fields are extinguished. The cloud is imaged by 
pulsing the MOT beams for 200 and recording the re- 
sulting fluorescence image (in the xz plane) using a CCD 
camera. Fig. [2] shows the measured temperature T as 
a function of the detuning of the MOT beams, S. The 
error bars represent uncertainty estimates from variation 
in the expansion velocity determined independently along 
the X and z directions. The upper solid curve shows the 
prediction of simple Doppler theory for a two-level atom 
[17| , while the lower curve is a fit to a characteristic sub- 
Doppler scaling law [H, [l^ 

where / is the total intensity at the MOT position and 
Is = 72 mW/cm^ is the saturation intensity. The result- 
ing fit gives C„+„- — 0.38 ± 0.02 and a predicted mini- 
mum temperature To — (77±10) /iK in the zero-intensity 



FIG. 2: (Color online) MOT temperature vs. detuning for 
I = 14.3 mW/cm^. The upper curve is expected for Doppler 
cooling a two-level atom, while the lower curve is a fit to 
Eq. IT}. The atom number ranged from a minimum value 
of approximately 10* at 5/27r = —5 MHz to a maximum 
of 1.3 X 10^ at S/2tt = -30 MHz (-O.SF). (Inset) MOT 
temperature vs. trapping beam intensity at a fixed detuning 
5 — — 0.7r. The shaded gray area is a fit to Eq. ([T]) together 
with the spread resulting from uncertainty in S. Error bars 
are described in the text. 

limit. We also varied the intensity at fixed detuning 
S = -0.7r (inset of Fig. [21), to find C^+^- = 0.15 ± 0.08 
and Tq — (125 ± 21)/iK. The difference in fit parame- 
ters between the data sets indicates that the tempera- 
ture may have a more complicated dependence on in- 
tensity than the simple scaling law of Eq. ([1]). This is 
not too surprising, since in our trap, intensity-dependent 
rates of optical pumping into magnetically trapped dark 
states may strongly affect the MOT temperature (we 
show below that the magnetically trapped atom frac- 
tion exhibits anomalous temperature behavior). Never- 
theless, our results show qualitatively and quantitatively 
that strong, efficient sub-Doppler cooling occurs in the 
MOT. Wc emphasize that this is a single-stage process 
(see also Ref. (3), where sub-Doppler coohng occurs in 
the MOT without an additional "optical molasses" pe- 
riod of polarization-gradient cooling in a nulled magnetic 
field. 

As discussed in Refs. [l^, [2l|, cr+cr^ polarization- 
gradient cooling typically breaks down in a longitudi- 
nal magnetic field because the "locking velocity," where 
atomic-motion induced frequency shifts are balanced by 
Zeeman shifts, is generally different for sub-Doppler and 
Doppler cooling mechanisms. This difference arises be- 
cause Doppler mechanisms rely on absorption and spon- 
taneous emission, so Zeeman shifts of both the ground 
and excited states affect the cooling efficiency; on the 
other hand, sub-Doppler mechanisms rely on stimulated 
processes between ground state sub-levels, and thus de- 
pend primarily on Zeeman shifts in the ground state. Be- 
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FIG. 3: (Color online) Semiclassical force F{p) (in units 
of hkr), momentum diffusion coefficient D{p) (in units of 
h^k^r), and equilibrium probability distribution W{p) (arbi- 
trary units) for ^^'^Rb (dashed lines, red) and ^^''Er (solid lines, 
blue) in a a'^a~ optical molasses in a longitudinal magnetic 
field B = 1 mT. For both cases, the saturation parameter 
s = 2 and detuning S = — 2.5r. 



cause the Lande g factors of the ground and excited states 
of a cooling transition are generally unequal, there is a 
critical magnetic field value at which the two mechanisms 
"unlock," with the larger capture range of the Doppler 
cooling force dominating beyond this point [l^, [2l| . For 
most laser-cooled atoms, this critical magnetic field value 
is around 0.1 mT. In Xe, for example, the difference 
in g factors between the ground and excited state is 
^Qeg = —0.17 and polarization-gradient cooling was ob- 
served to fail at 0.2 mT (2^. On the other hand, the 
g factor difference for the 401 nm transition in Er is 
Ageg = —0.004 [l3|, suggesting that cr+fT^ polarization- 
gradient cooling should remain effective in longitudinal 
magnetic fields as large as 10 mT. 

To investigate the possible effect of g factor degener- 
acy on Er sub-Doppler cooling in a magnetic field, we 
calculated the semiclassical force F{p), momentum dif- 
fusion coefficient D{p), and steady-state momentum dis- 
tribution W{p) for an atom with momentum hkp in a 
one-dimensional a'^a^ optical molasses in a longitudinal 
magnetic field B. Our calculations follow the method of 
Ref . t20j , in which p is treated as a classical variable and 
a Fokker-Planck equation is derived by expanding the 
optical Bloch equations to second order in hk. Results 
are shown in Fig. [3] for the 401 nm line of ^^^Er and for 
the D2 line of ^^Rb {Ageg = —0.17) in a magnetic field 
B — I mT. Because of the complicated interplay between 
optical coherences, Zeeman coherences, and atomic mo- 
mentum, the distributions W{p) are clearly not Gaus- 
sian. As a result, a temperature does not exist in the 
strict sense, but we may still take the variances of W{p) 
as a measure of the momentum spread. These variances 
correspond to center-of-mass "temperatures" of 25 /iK 
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FIG. 4: (Color online) CCD image of magnetically trapped 
atoms after 0.3 s in the magnetic trap. The marginal dis- 
tributions and fits to Px{x) and Pz{z) are shown at bottom 
and left, respectively. For this image, the fit parameters are 
X = 150 fim, z = 105 fim, and g = 0.23, corresponding to 
an effective magnetic moment p, = 3.9 yUs and temperature 
T = (44 ± 7) fiK. 



and 3.6 mK for ^^^Er and ^^Rb, respectively. The same 
calculation gives respective temperatures of 25 /iK and 
17 /iK at zero magnetic field B = Q. For this simplified 
one-dimensional case, the (/-factor difference results in a 
temperature increase of more than 2 orders of magnitude 
for *^Rb, with no change in temperature for ^^^Er, at a 
magnetic field value of i? = 1 mT (10 G). Note that this 
value of the magnetic field is attained in our setup at 
a distance of a few millimeters from the MOT position, 
within the beam waist of our trapping lasers. The near 
degeneracy of Er factors may contribute to capture and 
loading processes in our MOT and should also have in- 
teresting applications for laser cooling Er in moderately 
large magnetic fields. 



III. MAGNETIC TRAPPING 

We now consider the fraction of ground-state atoms 
that remain magnetically trapped after the MOT beams 
have been turned off. Assuming these atoms can be de- 
scribed by an effective magnetic moment /l > 0, then the 
potential energy for a trapped atom is given by 



U{x,y,z) = fiJdB'^x'^ 



dB^v^ 



Mgz 



where M is the atomic mass, and g is the gravita- 
tional acceleration. In equilibrium at temperature T, 
the marginal distributions of atoms along the x and z 
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FIG. 5; (Color online) Magnetic trap temperature vs. MOT 
temperature. The MOT temperature was varied by chang- 
ing the detuning S. Blue circles and red squares represent 
data from separate experimental runs at different Zeeman 
slower and MOT beam intensities. Simple thermodynamic 
arguments predict that all points should lie above the dashed 
hue where Tb > Tm /3. 



FIG. 6: (Color online) Evolution of the magnetically trapped 
atom number. After a time At « 0.2 s, all of the metastable 
atoms have returned to the ground state and the trap decays 
exponentially with a time constant of (270 ± 13) ms (solid 
line) . (Inset) The temperature of the trapped atoms increases 
over time with an error- weighted heating rate of (14 ±4) /iK/s 
(solid line). 



directions are found to be 



p,(z) = exp ( -2M „ 25^) (1 + 2M 
z z I \ z 



where Nx and Nz are normalization constants. The 
length scales and effects of gravity are characterized by 
g = Mg/{pbdBz) and x = 2kBT/{p,dBx) with a similar 
expression for z. Such a trap is stable against gravity only 
if ^ < 1. For g = 0, X and z are the standard deviations of 
the cloud, while for g > 0, the cloud stretches along grav- 
ity. We capture images of magnetically trapped atomic 
clouds by loading the MOT for 2 s then extinguishing the 
loading and trapping lasers, leaving all magnetic fields 
unchanged. After a variable hold time At, we pulse the 
MOT beams and capture the resulting fluorescence im- 
age. Such an image is shown in Fig. U] By analyzing 
these images, we can determine T, /2, and atom number 
as the magnetic trap evolves in time [l^] . 

If the magnetic trap loading efficiency is independent 
of an atom's energy, then thermodynamic arguments [2^ 
predict that the magnetically trapped atom temperature 
Tb should be related to the MOT temperature Tm by 
Tb > Tm /i (equality holds when the spatial extent of the 
MOT is much smaller than that of the magnetic trap). 
However, our magnetic trap temperatures systematically 
violate this bound as shown in Fig. [S] where we mea- 
sure temperatures below 25 /iK. Similar anomalously low 
temperatures have been observed in MOT-loaded, mag- 
netically trapped Ca [l^l and Sr [11] while unexpect edly 
high temperatures were reported in Cr [2^. In Ref. [25| |. 



a selection effect during trap loading was suggested to 
explain the low temperatures. We suspect a similar phe- 
nomenon is at work here, in which spatially-dependent 
optical pumping in the MOT tends to correlate an atom's 
spin state with its energy; the ejection of atoms in un- 
trapped magnetic states (mj < 0) may then reduce the 
average energy of the sample when the MOT is turned off. 
This hypothesis cannot be quantitatively evaluated with- 
out studying the loading efficiency of the magnetic trap, 
which is difficult to assess due to the presence of trapped, 
dark-state atoms that cannot be probed directly. 

An example of the time evolution of the magnetic trap 
population and temperature is shown in Fig. [5) Before 
Ai w 0.2 s, the trap population appears to increase. This 
is due to the relaxation of metastable atoms back to the 
ground electronic state where they are resonant with the 
probe light (MOT beams). After this time, we see single- 
exponential decay, most likely arising from collisions with 
the background gas, with a pressure around 2. 3x10^^ Pa 
[28| . This claim is further supported by a strong correla- 
tion between background gas pressure and trap lifetime 
and the similar lifetimes for ^^^Er and -'^^^Er (data not 
shown). We do not suspect that the trap decay rate 
is due to density-dependent coUisional effects within the 
trap, because the number densities in the present mea- 
surements are not particularly large [we observe peak 
densities of n = (2.1±0.3) x 10^ cm^"^] and because we see 
no evidence of non-exponential decay. On the other hand, 
we do not know the origin of the observed heating rate, 
which may arise from dipolar relaxation processes or col- 
lisional heating mechanisms. Improvements to the vac- 
uum system should eliminate our dominant loss channel 
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and allow us to study other relaxation processes. Finally, 
our data indicate a phase-space density p « 2.5 x lO^'*. It 
should be straightforward to increase this by at least an 
order of magnitude by increasing the atomic beam flux. 

IV. CONCLUSIONS 

In summary, we have presented a number of favorable 
cooling and trapping properties for atomic erbium. We 
showed that efficient sub-Doppler cooling occurs in our 
MOT without the need for an additional optical molasses 
cooling stage. Furthermore, sub-Doppler cooling should 
remain effective in large magnetic fields due to the near 
degeneracy of Lande g factors. This phenomenon may 
have interesting applications for laser cooling in large 
magnetic gradient traps, and might be useful for con- 
structing atomic waveguides along magnetic field lines. 



We also studied the thermodynamics of a continuously- 
loaded magnetic trap, where temperatures reach unex- 
pectedly low values. The eflticient cooling and trapping 
methods discussed here should be straightforwardly use- 
ful as a starting point for phase-space compression to- 
wards quantum degeneracy. However, magnetic and col- 
lisional relaxation processes must be studied before such 
a procedure can be considered. 
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